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We investigate separability and entanglement of mixed 
states in C 2 ®C 2 ®C N three party quantum systems. We show 
that all states with positive partial transposes that have rank 
< TV are separable. For the 3 qubit case (N = 2) we prove 
that all states p that have positive partial transposes and 
rank 3 are separable. We provide also constructive separa- 
bility checks for the states p that have the sum of the rank 
of p and the ranks of partial transposes with respect to all 
subsystems smaller than 15N-1. 

03.67.Hk, 03.65.Bz, 03.67.-a, 89.70.+C 



I. INTRODUCTION 

In the recent years it became clear that entanglement 
is one of the most important ingredients of the quantum 
information processing. While in the early age of quan- 
tum mechanics, entanglement was associated with "para- 
doxes" of quantum mechanics in the last decade 
of the last century it has been discovered that entan- 
glement plays an essential role in fundamental applica- 
tions of quantum mechanics to information processing 
(cf. |^-||). While the characterization of separable and 
entangled pure states of bipartite systems is quite well 
understood (cf. Q), it is not the case for mixed states. 
In the last four years, however, a lot of progress has been 
achieved in our understanding of the separability and en- 
tanglement problem for bipartite systems (cf. JtJ]). The 
first major step was the proper definition of separable 
and entangled states formulated by Werner || . The next 
milestone was the discovery by Peres || of the fact that 
all separable states are must necessarily have a positive 
partial transpose |ic| ]. Soon after Horodeckis |llj have 
shown the Peres criterium provides also a sufficient con- 
dition for separability in two qubit (2 x 2) and one qubit 
one qutrit (2 x 3) systems. Subsequently, P. Horodecki 
|jl2f has constructed the first examples of the, so called, 
bound entangled states, i.e. the first examples of the en- 
tangled states with positive partial transpose (PPT ES). 
This discovery has stimulated great interest in the stud- 
ies of properties of PPT ES. Some of the most important 
results, in particular coming from the Horodecki family, 
IBM group, and Innsbruck-Hannover collaboration are 
described in Refs. 0. 

More recently, considerable interest has been devoted 
to multiparty entanglement fl3[| . The first papers on 3 
qubit states led to the discovery of the, so called, GHZ 
states H|, which are particularly suited to study the 



break down of the Bell like inequalities in quantum me- 
chanics Q . Three party (and multiparty, in general) en- 
tanglement of the GHZ type can allow for interesting ap- 
plications, such as for instance quantum secret sharing 
JPoj ], and many experimental groups have recently tried 
to generate such states Q . 

Theoretical studies of the structure of multiparty en- 
tangled states has just started |l3| . First of all, pure state 
entanglement has been investigated. An important direc- 
tion of research was here initiated by Ref. flTf . In this 
paper Linden and Popescu have studied whether a given 
quantum state can be transformed into another one using 
local unitary (or at least non-unitary invertiblc) transfor- 
mations. Such a geometric approach calls for studies of 
invariants of local unitary and local non-unitary invert- 
ible transformations, and leads elegantly and naturally 
to the concepts of Schmidt coefficients ||, and Schmidt 
number |l8| for pure states in bipartite systems. 

This approach and concepts can be generalized to the 
case of 3 qubit systems and in general for mixed states, 
but it is by no means an easy task. In particular, as 
pointed out in Ref. in both cases one expects vari- 
ous, locally not equivalent kinds of entanglement to arise. 
Very recently the concept of Schmidt coefficients (i.e. in- 
variants of the local unitary transformations has been 
formulated for 3 qubit systems The other ap- 

proach (based on the investigations of local non-unitary 
invertible operations) has been followed by the Innsbruck 



group |22|. Dur et al. were able to show that there are 
essentially 3 types of entanglement of pure states: bi- 
partite entanglement, PF-class entanglement, and GHZ- 
class entanglement. The are ways of characterizing the 
3-qubit entanglement with the help of a, so called, tangle 
p3| , and other local invariants |2^,^l[. Numerous stud- 
ies of various types of multiparty entanglement of pure 
states and various interesting examples of it have been 
conducted in the recent years (24j| . 

At this point it is worth mentioning that while for the 
pure states in bipartite systems it is possible to quantify, 
or better to say to characterize the entanglement in the 
canonical way , this is not necessarily the case for the 
mixed states. The studies of entanglement and separa- 
bility in mixed states of a three party system has just be- 
gun. Among the recent results it is worth mentioning the 
demonstration of separability of the states that differ not 
much from fully chaotic state |2^] , the construction of en- 
tangled states that have all partial transposes positively 
defined, employing the concept of unextendible product 
basis pq|, the classification of multi-qubit states based 
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on the separability and distillability properties of certain 
partitions, pTj], the generalization of the concept of mean 
Schmidt number to the case of multiparty systems p8| , 
a formulation of the necessary and sufficient conditions 
for separability in term of linear maps ]29| , studies of the 
properties of relative entropy in multiparty systems |30| , 
and studies of states symmetric with respect to trilateral 
unitary rotations |3l| . We have presented recently a 
classification of mixed states for 3 qubit systems into the 
separable class, the bipartite, the W-, and the GHZ- 
classes of states. Following the Refs. |33|] we constructed 
canonical form of entanglement witnesses for each class, 
and discussed their optimization. 

In this paper we consider entanglement and separabil- 
ity of mixed states in C 2 ® C 2 (g> C N three party quantum 
systems. Such systems are of practical interests since i) 
for N — 2 they reduce to the intensively studied 3 qubit 
systems; ii) for N large, they can be used to describe 
two qubits interacting via a "bus" mode; this is how the 
quantum gates can be realized in the quantum computer 
model based on cold trapped ions . 

This paper generalizes the results obtained by us ear- 
lier for the case 2 x N |56| and M x N |37| systems. 
We use here the same mathematical tools that have been 
developed in our earlier work p8|j3S| , i.e. the method 
of subtracting from a given state p projectors on prod- 
uct states keeping the remainder, as well as its partial 
transpose (-s) positively definite 

The paper is organized as follows. In section 2 we 
demonstrate that all states in C 2 £g>C 2 ®C systems with 
positive partial transposes that have rank < N are sep- 
arable. This section is divided into 3 subsections, and 
the main result is presented in the last subsection. In 
the first subsection we present the canonical form of the 
investigated states; in the second one we prove an impor- 
tant Lemma that states that for the 3 qubit case (N = 2) 
all states p that have positive partial transposes and rank 
3 are separable. In the section 3 we discuss constructive 
criteria and separability checks for the states p that have 
the sum of the rank of p and the ranks of partial trans- 
poses with respect to all subsystems smaller equal than 
15./V — 1. We discuss here the concept of the edge states, 
i.e. those from which no projector on a product state 
can be subtracted without loosing either the positivity, 
or the PPT property. We discuss here also the methods of 
constructing the, so called, entanglement witnesses, and 
their canonical form. In section 4 we specify the previous 
results for the case N = 2, and we provide constructive 
separability checks for the states p that have the sum of 
the rank of p and the ranks of partial transposes with 
respect to all subsystems smaller equal than 29. 

In this paper we denote by R{p), K(p), r{p) and k(p) 
the range, the kernel, the rank, the dimension of the ker- 
nel of p, respectively. Also, |e) will denote a vector or- 
thogonal to \e). The symbol diagfci, <72, • • •] denotes a 
matrix with diagonal blocks o~i, o~%, 



II. PPT STATES OF RANK N IN C 2 ® C 2 ® C N 
SYSTEMS 

A. Generic form of the rank TV PPT states 

In this section we will derive the canonical form of the 
separable states in C 2 ® C 2 <g> C N with r(p) = N. The 
canonical form will allow for an explicit decomposition 
of a given state in terms of convex sum of projectors on 
product vectors. In the following the three parties will 
be called Alice, Bob and Charlie. We begin with the 
following Lemma: 

Lemma 1 : Every PPT state p in C 2 ® C 2 (g> C N 
with r(p) — N , such that in some local basis (\0a), 
\1 A ) for Alice, |0 B >, |1 B ) for Bob, \0 C ) . . . , \N - l c ) 
for Charlie) without loosing the generality we have 
r{{i-A, 1b\p\Ia, 1.b)) = N, can be transformed using a 
reversible local operation to the following canonical form: 



P = 



I B^C^CB B^C^C B^C^B B^C^ \ 
C^CB C^C C'^B Ct 
B^CB B^C B^B St 

V CB C B 1 / 

/ Btc* \ 

V 1 / 



(CB C B l)y/D, 



D 



(1) 



where [B, B^} = [C, Ct] = [C, B] = [C, fit] = o and D = 
B,C and D are operators acting in the Charlie's 
space. 

Proof: The state p can be always written in the consid- 
ered basis as: 



/ E 1 



El 

4 
\4 



E 5 
E-> 



E. 



Eg Ey \ 

^8 Eg 

E3 E10 



E l E\o E 4 



where E's are N x A^-matrices, and r(£vt) = N. After 
the projection p = (1a\p\1a) we obtain the reduced state 



P = 



E, 
E 



10 



E\Q 

E 4 



After performing a reversible local non-unitary " filtering" 
— A== on Charlie's side the matrix p can be written as: 

A fit 
B 1 

This matrix is obviously positive, i.e. can be represented 
as §7| p = £ + diag[A, 0], where A = A - B^B, 



E = 



( B^B St 
I B 1 



2 



The matrix p mufi has the rank N . We observe that 
E has also the range N, and possesses N vectors in its 
kernel \<f> f ) = - 12)51/). We will show that A = 0. 

Using the fact that p > 0, we observe that A > 0. 
But, since r(p) = r(E), the ranges of the the matri- 
ces must fulfill R(p) = i?(E) D J?(diag[A, 0]), so that 
the corresponding kernels fulfill K (diag[A, 0]) D 
The kernel K(E) is spanned by the vectors of the form 
\(j)f) = 1 — \2)B\f), where |/) is arbitrary, for which 
(0/|diag[A, O]|0/) = must hold also. This means, how- 
ever, that A|/) = for all |/), and thus A = 0. 

The fact that B is a normal operator follows from the 
fact that p tA must be positively definite. This condition 
implies that BB^—B^B > 0. The latter positive operator 
has, however, the trace zero, and must therefore vanish, 
i.e. [B,B*] = 0. 

Similarly, if wc consider the projection (1b\p\^b), for 
the same reasons as above we conclude that the resulting 
matrix 

/ etc* ct 

P={ C 1 

with [C, Ct] = 0. Summarizing, after performing a local 
filtering operation —= we can bring the matrix p to the 
form: 



/ E l 

4 
4 
K4 



c^c 

4 
ct 



E 6 e 7 \ 
E 8 Ct 
B^B B 
B 1 J 



Now, the matrix p possesses as kernel vectors |10)|/) — 
and |01)| fl ) - |H)C| 5 ) for all \f),\g) from' the 
Charlie's space. This implies that we must have E 8 — 
C^B, E 6 = E 7 B and E 5 = E 7 C. The matrix p has thus 
the form: 

/ E\ E 7 C E 7 B E 7 \ 

&e\ c^c &c ct 

b*e\ B^C B^B fit 



E] C B 1 / 



In the next step we consider its partial transpose with 
respect to Alice given by: 



4a 



/ E x E 7 C B^E\ B*C\ 

&e\ c^c e 



E 7 B 
V C*B 



E 7 

ct 



t 

7 

B*B 
B 



Ct 

st 
1 



Since partial transpose with respect to Alice is positive 
and does not change (I^IpIIa), the vectors |10)|/) — 
|11).B|/) should remain in the kernel. This implies the 
equality E 7 — B'C' , and the following form of p: 

( E x B^C^C B^C^B £tct ^ 

&CB CtC &B Ct 

B^CB B^C B^B St 

\ CB C B 1 J 



P 



The above form can be rewritten as 

f B^C^CB B^C^C B^C^B B^C^ \ 
C^CB C^C C^B Ct 
B^CB B^C B^B St 
CB C B 1 



+diag[A, 0,0,0], 



where A = Ei—B^C^CB. Using the short hand notation 
we get 



P = 



( Btct \ 
Ct 

1 



CB C B 1 ) +diag[A, 0,0,0]. 



The first term in p is PPT and has the following 3iV 
vectors in the kernel: 

|V) = |00)|/) + |ll)CB|/) 
\<j>) = |01)|. 9 ) + |11)C| 3 ) 
| X > = |10>|ft) + |ll)S|/»>, 

for arbitrary \f),\g) and \h). Similarly as above, this 
means that as in the case of A, the matrix A must vanish. 
This provides us with the final form of p: 



( B^C^CB B^C^C B^C^B B*& \ 

C^CB C^C C^B Ct 

StC B^B St 
C B 1 



B^CB 
CB 

( \ 

ct 
st 

1 



CB C B 1 



(2) 



(3) 



It remains only to prove the commutation relations 
[B, C] = [B, Ct] = 0. This follows from the positivity 
of all partial transposes of p. In particular, p tA is: 



p tA = 



C 

st 
1 



( &B Ct B 1 ) , 



(4) 



which is obviously positive definite. 
In contrast, p tB can be written as: 

/ B^t OB CtCB Bi&B &B\ 

fit CtC CtC BtCt Ct 

B^CB CB B^B B 
B^C C £t i 



(5) 



Because of its positivity, the matrix p tB must possess 
the kernel vector |01)|g) — |ll)C|.g), which implies that 
[C, B] = 0. The matrix p tB can be then written as: 



ct 

B 
1 



BtC C fit i ) 



(6) 
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which implies automatically the positivity. It remains 
finally to consider p tAB . The latter can be written as: 



/ B^&CB C^CB B^CB CB \ 

B^B C'^B B^B B 
\ B^C^ C+ fit i j 



(7) 



From the positivity of p tAB follows that |10> -\U)B^\f) 
is a kernel vector, so that [B\ C] = must hold. This in 
turn allows to write p tAB as: 



(CB\ 

C 
B 

V 1 / 



( Btct C+ fit i 



(8) 



Again, this form assures positive definiteness, and con- 
cludes the proof of the Lemma. □ 
Now, we are in the position to prove: 

Lemma 2 : A PPT-state p in C 2 ® C 2 <g> C N , whose 
rank r(p) = N , and for which there exists a product basis 
\&A,fB), such that r((e A , f B \p\e A , f B )) = N, is separa- 
ble. 

Proof: The state p can be written according to the 
Lemma ([l]) as 



/ Bttft \ 

V i J 



CB C B 1 



(9) 



Since all operators commute, they have to have common 
eigenvectors |/„), with eigenvalues b n ,c n , respectively, 
and 



(fn\p\fn) 



( K< \ 

K 

V i J 



Cnbn 



1) 



(10) 



This is, however, a product vector in Alice's and Bob's 
spaces. We can thus write p as p = J2n=i IV'nXV'nl ® 
\4>n)(4>n\ 8> \fn)(fn\- Because the local transformations 
used above were reversible, we can now apply their in- 
verses and obtain a decomposition of the initial state p 
in a sum of projectors onto product vectors. This proves 
separability of p, and the Lemma. □ 

From the Lemma 1 and 2 we conclude that in order to 
prove that PPT states p supported on C 2 ®C 2 ®C N with 
r(p) = N are separable, it is enough to show that one can 
find a product basis such that r{{e A , JbIpIs-a, /b)) = 
We will accomplish the proof in another way. Instead, 
we will prove the separability directly, and the desired 
canonical form of p will be a consequence of that. To this 
aim we will use the results of Ref . |36| , and the following 
theorem from Ref. |B7fl. 



Theorem 1 : For all PPT states p that are supported 
on M x N -space (M < N), and that have rank N, there 
exists a product basis such that r((lA\p\^A)) — N and p 
is separable and has the form 



N 



i=l 



,bi}(ei,bi 



(11) 



where are linearly independent. Additionally, the 
above decomposition is unique. 

This theorem can be used to prove the following 
Lemma: 

Lemma 3 : Any PPT state p supported on C 2 ®C 2 ®C N 
with N > A, for which r(p) — N , is separable, and obeys 
assumptions of Lemma 1. 

Proof: A C 2 ® C 2 <8> C^-system can be regarded as a 
C 4 ® C^-system. From the theorem [j] we obtain that: 

N 

P = \i>AB 1 XV'AB! | ® |Ci)(Ci | + W-ABi ' CiY&AB, ,Ci\. 

(12) 

Note, however, that we can find now a vector |C) in Char- 
lie's space, so that (C\p\C) ~ \ip ab^iP abA- Because the 
state p has a PPT property with respect to all partitions, 
li'AB^ii/JABi | must be PPT with respect to Alice or Bob, 
i.e. it must be a product state. This observation con- 
cerns all projectors that enter the convex sum (|l^), so 
that we conclude that p is separable. From (12) it fol- 
lows directly that p can be projected onto |1^, 1 B )> so 
that r((l A ,l B \p\l A ,l B ) =N.O ' 

B. Separability of states with rank < 3 in C 2 ® C 2 ® C 2 

Now we have to consider the cases N = 2,3. The 
following Corollary and Lemma deal with the case N = 2: 

Corollary 1 : Any PPT state p supported on C 2 ®C 2 ®C 2 
and such that r{p) = 2, has a product vector \e,f,g) in 
its kernel. 

Proof: The vector |e, /, g) belongs to the kernel iff it is 
orthogonal to the range, i.e. iff it is orthogonal to the 
two vectors {|^i), IV^)}, which span the range of p. We 
can choose |e) arbitrary and set |/) = |0) +a|l), so that 
we obtain two equations: 



((^|e,0)+a(^|e,l))| 5 ) =0 



(13) 



We treat these equations as linear homogeneous equa- 
tions for \g); they have nontrivial solutions if the corre- 
sponding determinant of a 2 x 2 matrix vanishes. This 
gives a quadratic equation for a, which has always at 
least one solution. □ 
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Lemma 4 : Any PPT state p, supported on C 2 ®C 2 ® 
C 2 , and such that r(p) — 2, is separable (compare 

Proof: If | eA, /s, gc) is in the kernel of p, then PPT 
property implies that also p tA \e* A , fs, gc) = 0. We 
obtain then that (e* A \p tA \e* A , f B , gc) = 0, where \e A ) 
is orthogonal to |e^). This equation is equivalent to 
(e A \p\e A , f B ,gc) = 0. This, however, means that 
p\eA,fB,gc) = \e-A,tpBc), where \ipBc) a vector in Bob's 
and Charlie's space. Now, according to the Lemma 2 
of Ref. |36| which deals with C 2 ® C^-systems, we can 
subtract the projector |e^, ipBc)(^A, ^bc \ from p, so that 

P = P~ (&a,1>bc\p-H&a,1>bc) 1 Ba ' ^Bc)(e A , V'bc I is positive, 
has rank 1, i.e. is a projector. Since it has the PPT prop- 
erty with respect to Alice's system, it must be separable 
with respect to A — BC partition. In general, we can thus 
write p = A|eAXe^| ® \'tpBc){'tpBc\ +^A,i>BC^fiA, ^sol- 
Projecting onto \ca) we get (e A \p\e A ) oc \tJjbc)(^bc\- 
Since p has the PPT property with respect to all par- 
titions, the projector IV'bcXV'bcI must project onto a 
product vector. The same can be of course said about 
I^'bcXV'bcIj since the projection onto |e^X^U| gives 
(e^|p|e A ) oc \iPbc){iPbc\, which implies that \ip B c){4>Bc\ 
is a product state and concludes the proof. □ 

Now we have still to prove the case N = 3. Before 
we do that, however, we need one more Corollary and 
Lemma concerning the case N = 2: 

Corollary 2 : Any PPT state p, supported on C 2 ® 
C 2 £§> C 2 , such that r(p) — 3, has a product vector |e, /, g) 
in the kernel. 

Proof: Let p be PPT-state in C 2 <g> C 2 <g> C 2 . It can be 
regarded as a C\ ® C|j C -state. According to Theorem 1 
of Ref. [[56| this state is supported on C A I S>C BC , and must 
have the form: 

3 

p = ^\e Ai )(e Ai \® lipBC^BCil (14) 

We take |e) orthogonal to \eA 3 ), and demand that |/, g) 
is orthogonal to \ipBd) and \ip B c 2 )- Setting \fs) = 
\0)b + ct\l)s, we obtain the following system of linear 
homogeneous equations for \g): 

{(il>BC<\0 B )+a(rl>BC t \lB))\g) =0 (15) 

for i = 1,2. These equations posses a nontrivial solution 
if the corresponding determinant of the 2x2 matrix van- 
ish. This lead to a quadratic equation for a, which has 
always a solution, and that proves the Corollary. □ 

The existence of product vectors in the kernel is used 
in the proof of the Lemma below. This Lemma provides 
one of the most important results of this paper: it implies 
that in C 2 (£)C 2 (£)C 2 systems there is no PPT entanglement 
of rank smaller than 4. 



Lemma 5 : Any PPT state p, supported on C 2 ®C 2 ®C 2 , 
such that r(p) = 3, is separable. 

Proof: Let \e A , f B , gc) belongs to K(p). From the con- 
dition p\eA, f B , gc) — for the product vector in the 
kernel, follows that: 

(e A \p\e A ,f B ,gc) = 0, 
(fB\p\e A ,fB,gc) = 0, 
(gc\p\eA,fB,gc) = 0. 

This means that 

p\e A) f B ,gc) = \e A )\ipBc), 
p\e A ,fB,gc) = \ fB)\i> A c), 
p\e A jB,gc) = \sic)\iPab)- 

We define 

p = p - \\e A )(e A | <g> | ipBCriipBC \, (16) 

where A = -r~ — ; — , Ii,- — ; — r (see Lemma 2 of Ref. 
HI). Now, p is a PPT state with respect to A - BC 
partition, i.e. p tA > 0; this state has the rank r(p) = 2. 
We rewrite p as: 

P = AiI/bX/bI <g) \iPac)(iPac\ + M\gc){gc\ ® \^ab)(^ab\- 

(17) 

We redefine now \ca) — |0) and \£a) — |1), i.e. change 
the basis in Alice's system, and represent the vectors 
\ipAc) and \4>ab) in the new basis as: 

|^c) = |0)|^) + |l)|V 2 ? ), (18) 

\^AB) = m l B )+\m%). as) 

In the matrix form p can be written as: 

/ I Ai|/b>(/b|® WcXV-fel A ( A 1 |/ B >(/ B |«.|^)<^| A \ 
-_ \ +A 2 |^><^|®|ff C Xsc| ) \ +*2\<l> 1 B ){<P 2 B \®\gc){gc\ J 

P ( a 1 |/bX/b|®I^XV'cI \ / AiI/bX/bI^I^cXV-cI \ 
\ \ +A2|0|)<0 B | ® \gc){gc\ ) \ +\2\4> 2 b){4> b \ ® \gc){gc\ J ) 

From the positivity of p and p tA follows that when a diag- 
onal block ( A il/sX/s| ® IV'cXV'cl ^ ortine-ou U 2 ib 2 ) 
onal block +M4>lm | g, |^ )( . c | J acting on |0 B V C ) 

vanishes, the same must be true for the off-diagonal block 
in left lower corner. Similarly, the same observation con- 

cerns the diagonal ( ^'{fj^lf ^ ){ f?\ Y the vec- 
\+^ B )^B\®\9c)(gc\ J' 

tor and the off-diagonal block in the right upper 

corner. This leads to the system of equations: 



(fB\ft B )^c\i>c) 


= 0, 


(20) 


(4>Wb)(9c\^c) 


-o, 


(21) 


(fB\$ 2 B )(^c) 


= 0, 


(22) 


(^ 2 b)(9c\^c) 


= 0. 


(23) 
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This system of equations implies that at least one of 
the projectors \4> ab)^ ab\ and \ipAc)(^Ac\ must be a 
product state. If it is, for instance, |V , ab)(V , ab| ) then 
\4> B ) = 10b) = I/b) and P becomes 

P = Ai|/b)</b| O IVacXV'acI 
+ A 2 |g A )(g j4 | <g> I/bX/bI ® IMsd 
+ Aje><e| ® I^bcXV'bcI 

= I/bK/bI ® (Ai|Vuc)<^ac| + A 2 |e^Xe A | <8> IscXsd) 



+ Aje><e| ® \4>bc){iI>bc\- 



(24) 



The operator tr is a PPT state of rank 2 in C 2 ® C 2 space 
of Alice and Charlie. From Peres-Horodecki criterium 
P.pr| follows that it is separable. The matrix p can thus 
be written as 

p = Ai|e A )M ® |/b)</b| (8 IscXffcl 
+ Aa^M 8> |/b)(/b| ® IscXffd 
+ A|e)(e|<8 I^bcXV'bcI. 

For the above proof Alice is in no way distinguished. We 
can also write 

p = p- X\e A )(e A \ <8> |/bX/b| <8> |ffcX.9c| 
= Ai|eAXe A | ® |/bX/b| ® IscXsd 
+ A|e A )(eA| <8> IV'bcXV'bcI, 

-, and p is a PPT 



where A = A 2 = - . _ , . 
state with respect to C—AB partition. The projection of 
p onto |e>t) gives (eA\p\e-A) ~ IV'bcXV'bcI- This means, 
however, that IV'bc) must be a product vector, and that 
concludes the proof. □ 



C. Separability of states of rank N in C 2 ® C 2 <g> 
systems 

Now we are in the position to prove the main theorem 
of this section. Before that we have to complete, however, 
the discussion of the case N = 3. To this aim we prove 
the following Lemma: 

Lemma 6 : Any PPT state p, supported on C 2 <g>C 2 <g>C 3 , 
such that r(p) — 3, is separable. 

Proof: We consider the system C 2 (g> C 2 (g> C 3 -System as 
a C AB €5 Cq system. According to the Theorem (|l|) three 
possibilities may occur: 

• The state is supported on C\ B <g> Cq. Then the 
density matrix must have a form 

p = Ki\e ABl ){e AB A ® 1/dX/cJ 
+ A 2 \e A B 2 )(eAB 2 \ ® \fc 2 )(fc 2 \ 
+ &3\eAB 3 )(eAB 3 \ <8> \fc 3 )(fc 3 \- 



Since the vectors \fc t ) are linearly independent, we 
can find a vector \C) in Charlie's system such that 
(C\p\C) ~ |e^BiX e ABil- Because the considered 
state has the PPT property with respect to all par- 
tition, the projected state |eABiX e ABi | is also PPT, 
and as such must be a product state. 

• The state is supported on C\ B ® Cq. The same 
method of projecting onto appropriately chosen 
vector in Charlie's space allows to prove the sep- 
arability. 

• The state is supported C^ s <X>Cp. That is, however, 
nothing else but a state in C 2 ®C 2 ®C 2 system with 
rank 3. Its separability follows from Lemma 5. 

This concludes the proof of the Lemma 6.D 

Now, all the above presented results can be brought 
together in a form of the following theorem: 

Theorem 2 : Every PPT state, supported on C 2 ® 
C 2 ® C , such that r(p) = N, is separable and has the 
canonical form of the Lemma [J, 

/ B^&CB B^C'^C B^C'^B S^t \ 
C^CB &C C^B Ct 
B^CB B^C B^B fit 
\ CB C B 1 / 

/ stct \ 
V 1 J 



CB C B 1 ) VD, 



D (25) 



(26) 



where B,C and D are operators acting in Charlie's space 
that fulfill [B, fit] = [C, C*] = [C, B] = [C, B*] = and 
D = Dt. 

In the next section we will study states in C 2 (g)C 2 (g)C N 
with low ranges, but > N. By looking at product vectors 
in the ranges of p and its partial transposes it is posiible 
to check separability for low rank matrices, similarly as 
in the case of bipartite systems in C M ® C N |pq,p7[. 



III. SEPARABILITY CHECKS AND CRITERIA 
FOR GENERIC LOW RANK STATES IN 
C 2 <g> C 2 ® C N SYSTEMS 

In this section we will study the PPT states p that 
posses a finite number of product vectors in their range 
\ei,fi,9i) ^C 2 ®C 2 ®C N such that \e Ai , /s 4 ,5c? 4 ) & R(p), 
le-XJn^gCi) € R{p tA ), \eA.,f Bi ,g Ci ) e R{p tB ) and 
\ e A>f B >9Ci) £R(p tAB ). We will show that this is gencr- 
ically the case when r(p) + r(p tA ) + r(p tB ) + r(p tAB ) < 
15 A — 1. Let us call the set of such vectors V[p] The 
search for the desired product vectors {\eA i: fBt, gd)} € 
C 2 ®C 2 ® C N is reduced to the problem of solving a sys- 
tem of multipolynomial equations [ ^6[ . When the num- 
ber of equations is equal to (bigger than) the number 



G 



of available unknown parameters, one expect the num- 
ber of solutions to be finite (zero). The states of low 
ranks fulfilling this property has been termed generic in 
Ref. |37|. In particular, the states for which the number 
of the desired vectors in any of the considered ranges is 
smaller that the corresponding rank, must be entangled. 
Particularly important are states that do not contain any 
product vector of the above described properties in the 
range. Such states are termed edge states, and play ma- 
jor role in characterization and classification of the PPT 
entangled states Q. 



A. Generic states 

Let \Ki), \K Ai ), \KBi) and \KABi) are linearly inde- 
pendent vector that span the kernels of p, p tA , p tB and 
p tAB , respectively, so that: 

K(p) = spsm{\Ki),i = 1, . . .,k(p)}, 
K{p tA ) = span{|A^), i = 1, . . . , k(p tA )}, 
K(p tB ) = span{|K Bi ), i = 1, . . . , k{p tB )}, 
K(p tAB )=sp&n{\K ABi ),i = l,...,k(p tAB )}. 

Choosing an orthonormal basis in Alice's and Bob's space 
we can write those vectors as: 



\Ki) 


= |oo)|4°M 


-|oi)|fc 01 ) + 


|10)|fcf) + 


in) 14) 


\K Ai ) 


= 100)1^°) H 


-101)14') + 


110)14°) + 


111)14, 


\K Bi ) 


= |oo)l4°H 


-101)14,) + 


|io)l4°) + 


111)14, 


\K ABi ) 


-|oo)|4° Bl ) 


+ 101)14%: 


> + |io)|4%) + in) 



A product vector in |e, /, g) € V[p] has the property that 
it and its partial complex conjugates have to be orthog- 
onal to the corresponding kernels, i.e.: 

(Ki\e A , fB,9c) = 0. 
(K Ai \e A ,f B ,g c ) = 0, 
(K Bi \e A ,f B ,g c ) = 0, 
(K ABl \e A ,f B ,g c ) = 0. (27) 

We expand now \e A , f B ,gc) in the local basis of Alice 
and Bob: 

\e A ,f B ,gc) = HO) + |1)) ® CS|0> + |1)) ® \g) 

= (a/3|00) + o|01) + P\10) + |11)) <8 \g). 

We observe that Eqs. ( p?| ) can be rewritten as : 

A(a,0;a*,F)\g) = O, (28) 

where A(a, j3; a*, (3*) is a (k(p) + k(p tA ) + k{p tB ) + 
k(p tAB )) x N matrix, which reads: 



A{a,/3; a*, (3*) 



a *f3(k 
cf3*(k 



:-:!! 
Ob 



+ affc? 1 ! + 3<fc! u + (*, 
+ ct*(k<* I + 3<fc™. I + <fc\ ] 



till 



11 1 



+ c*(fc 



01 



Eqs. ( |27| ) have a nontrivial solution with \e) 7^ 0,|/) =/= 
and \g) 7^ iff the rank of A is smaller than N. That 
implies that at most N — 1 rows of the matrix A are 
linearly independent. That means that (k(p) + k(p tA ) + 
k(p tB ) + k(p tAB ))-N+l minors of dimension N x N of 
the matrix A must vanish. 

Let us consider the marginal case, when k(p) + k(p tA ) + 
k(p tB ) + k(p tAB ) = 2 + (N - 1). In this case we combine 
the first N — 1 rows with the remaining two and obtain 
exactly two different minors, and thus two equations for 
complex a, j3, or more precisely four real equation for real 
and imaginary parts of a, (3. Such equations generically 
will have a finite number of solutions. The case when 
k(p) + k{p tA ) + k(p tB ) + k{p tAB ) > 2 + (N - 1), i.e. 

(r{p) + r(p tA ) + r{p tB ) + r(p tAB ) < 15N — 1 (29) 

means that we have more equations than parameters, and 
generically there will be no solution, or at least the num- 
ber of solutions will be even more limited than in the 
marginal case. The PPT states fulfilling the inequality 
(p9| ) are generically the edge states, provided their rank 
and/or the ranks of their partial transposes are greater 
that N, since otherwise the Theorem of the previous sec- 
tion would apply. Conversely, if 

(r(p) + r{p tA ) + r(p tB ) + r{p tAB ) > 157V -1, (30) 

then the matrix A has less equal than N rows, and one 
can always use the freedom of parameters to find a solu- 
tion, and subtract a projector onto a product vector from 
p keeping its positivity and PPT property intact. 

In the following we will concentrate ourselves on the 
case k(p) + k(p tA ) + k(p tB ) + k(p tAB ) >N+1, for which 
the number of solutions is expected to be finite. Such 
states will be called as in Ref. j}7j generic. For those 
states it is simple to check the separability, similarly as 
discussed in Ref. |^ , ^7| . The check is easy, because we 
know that if the considered state is separable, then it is 
represented as a convex sum of projectors on the vectors 
from the set V [p] , and the latter has a finite cardinality. 
We will discuss this in more detail for the case of C 2 £g> 
C 2 <X> C 2 systems. 



IV. SEPARABILITY CHECKS AND CRITERIA 
FOR GENERIC LOW RANK PPT STATES IN 
C 2 ®C 2 ® C 2 SYSTEMS 

As a special, but important example we consider the 
case of PPT states in C 2 ®C 2 ®C 2 states (3 qubit systems). 
We will use here the results of the previous sections. The 
3 qubit case is particularly interesting as a first step to- 
ward multiple entangled systems, providing a challenge 
for both the theory and experiment. 

Generically, if 



AB, 
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(r(p) + r{p tA ) + r{p tB ) + r(p tAB ) < 28, 



(31) 



E. The case r(p) = r(p tA ) = r(p tB ) = r{p tAB ) = 7 



then the set V[p] is empty and the state p is a PPT 
entangled edge state, provided all the ranks are greater 
3, since otherwise the Lemma 5 of the previous section 
applies. We discuss the different cases below 



A. The case r(p) = 2,3 

From the results of the previous section we know that 
such PPT states are separable. 



B. The case r(p) = 4 

The state of rank 4 in a 3 qubit system may be regarded 
a state in C 2 £g>C 4 of rank 4. From the Theorem [l] that this 
state is bipartite separable, and moreover has a unique 
decomposition into a sum of four projectors on product 
(biseparable) vectors in C\ ®C% C . From uniqueness, 
we gather that p is then separable iff the product vectors 
in this decomposition are completely separable, i.e. are 
product vectors in C 2 ®C 2 ®C 2 . Otherwise, the state is 
entangled, although biseparable. In fact it must be bisep- 
arable with respect to all partitions, i.e. also C 



'I' 



r 4 



and Cq ® C AB . Examples of such states are known, 
in particular those are the state constructed from unex- 
tendible product basis pq. 



C. The case r(p) = r(p tA ) = 5 

This case is also easy because first of all the bipartite 
separability with respect to the partition A — BC the 
has to be checked. As shown in Ref. J3(|. a PPT state in 
C\ ® C BC is (bipartite) separable, iff the set of bipartite 
product vectors Va-bc[p], corresponding to the partition 
A — BC is not empty. In the present case it must not 
only contain a bipartite product vector, but a tripartite 
product vector. If such vectors exist, generically there 
will be finite number of them, and at least 5 of them 
must belong to the set Va-b-c[p\- 



D. The case r(p) + r(p tA ) + r(p tB ) + r(p tAB ) < 28 

In this case we have more equations than available pa- 
rameters, and we expect that the set V[p] will be empty, 
whereas the state p will be an edge state. If this is not 
the case, we expect first of all that there is a finite num- 
ber of product vectors in V[p\. Thus, checking if p can 
be represented as a convex sum of projectors onto the el- 
ements of V[p] can be performed exactly using the same 
methods as discussed in Ref. ||36||. 



If V[p] is empty, this case describes an example of an 
edge state with maximal sum of ranks. Such an example 
has been constructed in Ref. |32|. Let us estimate how 
many elements can the set V[p] contain maximally. To 
this aim we write the matrix A: 



A(a,/3;a*,P*) 



°*' 3 < fc A 1 l + °*<4 1 1 l + ' j<fc Al l + <fc "l 



Let us denote a polynomial P of orders X and Y in vari- 
able z and z* by Px.y(z). Combining the first and the 
third row, and the second and the fourth row of A we 
obtain the minors of the form: 

a 2 P^ 1} (/?) + oP ( W } OS) + P ( W j 03) = 0. (32) 
(a'J^j^jOaj + o'flj^jOa) + R$ tl) (P) = 0. (33) 
The remaining four combinations of rows give us 

03) + «Q$o) (A + a *«(2!o) O 9 ) + Q(2?o) O 8 ) = ( 34 ) 

aa * Q(o?2) (/») + Q «(o! 2 ) 03) + a * Q\o,2) W + «(o?2) 03) = 0, (35) 
and two equations of the form 

<x<x*qV 1} 03) + aQg^ 09) + cfQWq 08) + Q$i) (/3) = 0. (36) 

Only 3 of the above equations are independent, but we 
have to our three complex conjugated equations_tq our 
disposal, and in particular the conjugate of Eq. (33), 



a 2 ii«* 03) + aH£J* 03) + «W* 09) = 



(37) 



A good strategy is to multiply Eq. ( j32| ) by R^l) (P) , and 



Eq. (|37j) by pi 1 ^, and subtract one from another in 
order to obtain 



a = T ( ( 2 1) 2) (/3)/T« ) (/3). 



(38) 



Inserting this solution into Eq. ( p2|) we obtain a polyno- 
mial of orders 5,5 in (3 and (3*. Another independent 
polynomial is obtained by complex conjugation. The 
variables (3 and j3* are then treated as independent ones, 
similarly in the Appendices of Ref. (36). According to 
the result presented there, a system of two polynomial 
equations of order X, Y with X < Y for two variables 
(3 and, say, (3 has at most 2 X Y solutions for (3. In the 
present case we expect thus that the number of solutions 
is < 160. Most of these solutions will have to be re- 
jected typically, since they do not fulfill the conditions 
Eqs. ©-©. 



F. The case r(p)+r(p tA ) +r{p tB ) +r(p tAB ) = 29 

This is a marginal case in which the number of equa- 
tions is equal to the number of parameters, so that gener- 
ically we have a finite number of product vectors in V[p], 



8 



and a possibility of performing the relatively straight- 
forward separability check. For example, if we consider 
r(p) = r(p tA ) — r(p tAB ) — 7 and r(p tB ) — 8. In this case 
only two minors are independent, and we have, for in- 
stance, to solve Eq. (|34|), one of the Eqs. ( |36j) , and their 
complex conjugates. By multiplying Eqs. (|34|) and j36| ) 
and its complex conjugates them by appropriate poly- 
nomials in /?,/?*, and subtracting one from another we 
obtain two linear equations for a, a* of the form 

aS$ A) 08) + a* S$ tl) 09) + sf A) (J3)=0. (39) 

and the complex conjugate of the above Eq. ([39]). This 
system of two linear equations can be solved so that we 
obtain 

« = r S)(/3)/rS)(/3)- (40) 

Inserting this solution into Eq. ( pf ) and obtain in this 
way a polynomial of order 9 in (3 and /3* . Another inde- 
pendent polynomial is obtained by complex conjugating 
Eq. (p^). The variables /3 and (3* are then treated as in- 
dependent ones, similarly as discussed in the Appendix 
of Ref . |3(| . According to Ref . we expect in this case 
maximally 2 9 x 9 = 4608 solutions for /3. 

G. Canonical form of non-decomposable 
entanglement witnesses 

For completeness it is worth mentioning that it is pos- 
sible to generalize the results of Ref. Q to case of 3 
qubit systems (and in general in tripartite systems). Let 
us remind the readers that an entanglement witness is 
a hermitian operator W, for which Tr(Wa) > for any 
separable state a, whereas Tr(Wp) < for some entan- 
gled state p. We say that W detects then p. A non- 
decomposable witness is a witness that detects a PPT 
entangled state. Using exactly the same arguments as in 
Ref. |33| one shows that a non-decomposable entangle- 
ment witness must have the canonical form 

W = P + Q tA + R tB + S tAB - el, (41) 

where 

e= inf (e,f,g\P + Q tA + R tB + S tAB \e, f,g), (42) 

|e,/,s) 

the operators P, Q, R, S are positively definite, R(P) = 
K(6), R(Q) = K(5 tA ), R(R) = K(5 tB ), R(S) = 
K(S tAB ), and 5 is an edge state, i.e. such state for which 
by definition the set V[5] is empty, which implies au- 
tomatically that e is strictly positive. According to the 
results of this section, in three qubit system, the state S is 
a generic state with r(5) + r(S tA ) + r(S tB ) + r(5 tAB ) < 28. 



V. CONCLUSIONS 

We have generalized previously obtained results for 
PPT state in C 2 (g>C N and C M ®C N system to PPT states 
in C 2 ®C 2 ® C N . We have developed a method of "local 
projections" together with the PPT property to prove 
separability of low rank states and to obtain separability 
criteria for low rank states. These methods together with 
methods developed in Refs. |^,^| provide very general 
mathematical tools to study separability and entangle- 
ment in multipartite systems. 

The main results of this paper are: 

• The proof that all states with positive partial trans- 
poses that have rank < TV are separable, and have 
a certain canonical form; 

• The proof that for the 3 qubit case (N — 2) all 
PPT states p that have rank 3 are separable; 

• The presentation of constructive separability 
checks for the states p that have the sum of the 
rank of p and the ranks of partial transposes with 
respect to all subsystems smaller than 15 N — 1. 

• The detailed discussion of the above mentioned 
constructive separability checks for the case N = 2; 

• Presentation of the canonical form of non- 
decomposable entanglement witnesses in 3 qubit 
systems. 

This work has been supported by the DFG (SFB 407 
and Schwerpunkt "Quanteninformationsverarbeitung" ), 
the ESF PESC Programm on Quantum Information, and 
European Union 1ST Program EQUIP. 
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